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Topological properties of QCD with two dynamical fermions 
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We investigate the topological susceptibility of the QCD vacuum with two flavours of dynamical staggered 
fermions on the lattice both at zero and finite temperature. At zero temperature we study the dependence of the 
signal on the fermion mass and at finite temperature we analyze the behaviour across the phase transition. 



1. INTRODUCTION 

The topological susceptibility is an important 
parameter of the QCD vacuum. In the continuum 
it is defined as 



(1) 



X = J d^x d^{0\T{K^{x)QiO)}\0) 
where K^{x) is the Chern current 



K^.{x) 
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(2) 



and Q{x) — dfj,K^{x) is the density of topological 
charge, 
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e'^''"^ F^(x) F;Ax) 



(3) 



Eq. (|^) uniquely defines the prescription for the 
singularity of the time ordered product when 

We present results on the topological suscepti- 
bility in QCD with Nf — 2 degenerate flavours. 

1.1. The Simulation 

We have simulated the theory on a 32'^ x 8 lat- 
tice with two flavours of staggered fermions with 
bare mass am = 0.0125. The R-type HMC algo- 
rithm has been utilized for the updating |^ . Each 
trajectory consisted in 60 steps of Ar — 0.005 
units of molecular dynamics time. We have per- 
formed the simulation at the following values of 



the lattice bare coupling /3 = 6/(7^ = 5.4, 5.5, 5.6 
and 5.7 . To be sure that the topology is well ther- 
malized we have checked the topological charge 
of our sample of configurations by cooling 
In Fig. |l| we show the distribution of topological 
charge from our configurations aX (i — 5.5 . The 
histogram displays a good sampling. To achieve 
the same ergodicity at larger values of j3 we had 
to increase the statistics. 

We have measured the Polyakov loop at zero 
and finite temperatures. With our bare fermion 
mass and lattice size the deconfining transition 
seems to occur at (3c = 5.54(2) 
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Figure 1. Distribution of topological charge after 
30 cooling steps at /3 = 5.5. 
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Table 1 

Topological susceptibility, T/Tc and a? a vs. j3. 







T/T, 


10-« x/MeV* 


5.4 


0.138(8) 


0.707 


5.16(1.66) 


5.5 


0.082(9) 


0.917 


6.21(1.74) 


5.6 


0.053(3) 


1.143 


2.27(0.72) 


5.7 




1.430 


1.72(1.00) 



1.2. The Operators and their Renormaliza- 
tions 

We have measured the topological charge den- 
sity on the lattice by using the operator Ql{x) M 



1 



Tr{U,,,{x)Up,{x)) , (4) 



and applying two smearing steps on it. The lat- 
tice topological susceptibility is defined by 
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(5) 



where Ql is the total topological charge and V 
the space-time volume. To match this definition 
to the one given in Eq.(|l]) we must subtract the 
singularity in the product of the two operators in 
Eq.(^). We make the subtraction by using the 
field-theoretical method M 



XL 



(6) 



and evaluating the multiplicative and additive 
renormalizations Z and M by use of the heating 
method Q . M contains mixings with gluonic and 
fermionic condensates. The lattice topological 
charge Ql mixes with fermionic operators dur- 
ing renormalization but it has been shown that 
the off-diagonal mixing is negligible . 

The heating method to evaluate Z and M 
yields a non-perturbative determination of these 
renormalization constants To calculate Z a 

local updating algorithm is applied on a config- 
uration containing a charge -1-1 classical instan- 
ton. These updatings, being local, thermalize the 
short distance fluctuations, responsible for the 
renormalization effects, and due to the slowing 
down leave large structures, like instantons, un- 
changed. The measurement of on such up- 
dated configurations yields Z -Q. As Q is known, 
one can extract Z. Notice that this procedure is 



equivalent to imposing the continuum value for 
the 1-instanton charge (in the MS scheme it is 
+1) and extracting the finite multiplicative renor- 
malization Z by evaluating a matrix element of 
Ql- 

The additive renormalization M is obtained in 
a similar way. We apply a few heating steps with 
a local updating algorithm on a zero-field con- 
figuration. Then the topological susceptibility is 
calculated. This provides the value of M, if no in- 
stantons have been created during the few updat- 
ing steps. This method agrees with the treatment 
of the singularity of Eq. (|l|) . 

As explained in [|| , cooling tests must be done 
to check that the background topological charge 
has not been changed during the local heating. 

1.3. Determination of the scale 

The lattice spacing was extracted by measuring 
the string tension at /3 = 5.4, 5.5 and 5.6 on a 
16^ lattice and assuming ^/a = 420 MeV. For 
f3 = 5.7 it was determined by an extrapolation 
using the 2-loop beta function. Wilson loops were 
evaluated by using smeared spatial links. The 
determinations of a^a are shown in Table |^. From 
the value at /3c (a = 0.123(6) fm) we infer the 
critical temperature to be Tc = 200 ± 10 ± 15 
MeV where the first error is our statistical one 
and the second comes from the indetermination 
in the result for Pc- 

2. RESULTS 

The topological susceptibility at zero temper- 
ature is shown in Fig. ^j. In the lower (upper) 
axis we show the values of (3 (bare fermion mass) . 
It displays a trend consistent with the theoretical 
expectation x oc m{^'ip). 

The values for the topological susceptibility at 
finite temperature are given in Table In Fig. | 
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we show the normahzed topological susceptibility 
x{T)/x{T = 0) as a function of the temperature. 
The value at zero temperature for (3 = 5.7 has 
been obtained by extrapolating the results shown 
in Fig. H In Fig. ^ the results for the quenched 
case Q are shown for comparison. The signal 
for X drops when crossing the transition temper- 
ature. Within the large errors, the drop seems 
to be as sharp as it was for the quenched case. 
Our results for Nf = 2 are in qualitative agree- 
ment with the ones shown in For Nf = 4 
the drop in the signal for x turns out to be much 
steeper [^21. 
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Figure 2. Topological susceptibility as a function 
of the bare fcrmion mass. 
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